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Abstract 

We propose a Lax equation for the non-linear sigma model which leads directly to the conserved 
local charges of the system. We show that the system has two infinite sets of such conserved 
charges following from the Lax equation, much like dispersionless systems. We show that the 
system has two Hamiltonian structures which are compatible so that it is truly a bi-Hamiltonian 
system. However, the two Hamiltonian structures act on the two distinct sets of charges to give the 
dynamical equations, which is quite distinct from the behavior in conventional integrable systems. 
We construct two recursion operators which connect the conserved charges within a given set as 
well as between the two sets. We show explicitly that the conserved charges are in involution 
with respect to either of the Hamiltonian structures thereby proving complete integrability of the 
system. Various other interesting features are also discussed. 
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1 Introduction: 



The non- linear sigma model has been studied from various points of view |ffl]. In the context of 
1 + 1 dimensional integrable systems, it is one of the most enigmatic models. The equations for the 
principal chiral model, for example, take the forms 

d^ a = 0, 

d^-d^-r^j'j^o, (i) 

where fx = 0,1 and belongs to a simple Lie algebra, with a taking values in the adjoint repre- 
sentation. Explicitly, the two equations can be written as 

d t J$ = QbJ? 

d t J? = d x J$ + f abc J*JZ (2) 

It is well known that the equations (|]) can be combined into a single zero curvature condition 
by introducing a one parameter family of connections that depend on a constant spectral parame- 
ter Q. The zero curvature representation is quite useful and leads directly to the infinite number of 
conserved non-local charges of the system ||] . While an infinity of conserved charges is a hallmark 
of integrable systems, these non-local charges are not in involution with respect to the Hamiltonian 
structure of the system Q. Rather, they satisfy a Yangian algebra p]. While classical integrabil- 
ity of the non-linear sigma model is known from alternate arguments |], 0] (see |§ for quantum 
integrability) , a conventional proof of integrability through the existence of an infinite number of 
conserved charges in involution is lacking so far. Furthermore, it is not known whether the system 
is bi-Hamiltonian, as is the case with most integrable systems || [Tcf| . 

There have been various attempts [11] to construct conserved local charges for the system (|l]), 



but a systematic study through a Lax representation is lacking so far. In this letter, we will provide 
a Lax equation for the non-linear sigma model equations @). Although the equations (JlJ) cannot 



be characterized as dispersionless, they share many features common to a dispersionless model [12 
In section 2, we discuss briefly the two component Riemann equation. This is a dispersionless 
system sharing some properties with the non-linear sigma model. In section 3, drawing from our 
experience in the Toda lattice |1C, 13], we present a Lax representation for (|l]). Even though the Lax 
equation for the system is not that of a dispersionless system, it immediately leads to two infinite 



sets of local conserved charges, a common feature of dispersionless systems [14, 15|. We point out 
various properties of these conserved charges that are quite useful in proving the integrability of 
the system. We also present the two Hamiltonian structures which are compatible and make the 
system bi-Hamiltonian. In section 4, we construct two recursion operators for the system which 
connect charges within any given set as well as between the two sets. We also show explicitly 
that the two sets of charges are in involution with respect to either of the Hamiltonian structures 
thereby proving the integrability of the system. In section 5, we summarize our results and discuss 
some other interesting aspects of this system. 
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2 Matrix Riemann equation: 



Before discussing the Lax equation for the non-linear sigma model, we will briefly discuss, in this 
section, the two component Riemann equation. We recall that the Riemann equation 

u t = uu x , (3) 

is an integrable dispersionless system [14, 16] which belongs to the class of equations of hydrody- 



namic type. Even though the non-linear sigma model is not a dispersionless system, it shares some 
nice properties with the two component Riemann equation. Let us define a 2 x 2 matrix 

U=( U V )=ul + va 1 , (4) 

where u, v are the two dynamical variables of the system and a\ is the familiar Pauli matrix (/ 
denotes the identity matrix.). Then, the generalized two component Riemann equation takes the 
form 

U t = UU X . (5) 

Let us next consider a matrix Lax function of the type 

L = p 2 I + U , (6) 

where p denotes the momentum in the phase space. Then, from the special structure of U in (]!]), 
it is clear that 

[a 1 ,U]=0, [a 1 ,L) = 0. (7) 
As a result, there is an ambiguity in taking the square root of the Lax function, namely, 

l\ =pl + hjp- 1 - l -U 2 p-* + --- . (8) 
Furthermore, a second square root has the form 

4 = *i4 • (9) 

We note that the square root (either one) of the Lax function commutes with cii, as does L. Both 



the square roots lead to consistent Lax equations 15] as 

dL (2n)H 



dt n (2n + l)!! 

dL _ (2n)!! 
d1~ ~ (2n + l)!! 



2n + l ~| 

(£i)>i 2 ,Lj , 

2n + l ") 

(^2)>! 2 ,L , (10) 
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where n = 0, 1, 2, 



and which lead to the hierarchy of Riemann equations of the form 



U t = U n U x , (11) 
U t = axU n U x , (12) 

respectively. In the second equation the effect of oi, on the right hand side, is to interchange u *-* v, 
and eq. (12) can be thought of as the "elastic medium" equations jl4|] associated with the hierarchy 
of Riemann equations (11). It is worth noting here that, although we have written the equations 
as non-standard equations, the same equations also arise from the standard Lax equations [17]. 

The conserved charges of the system can now be easily obtained from the Lax function. In fact, 
there are two infinite sets of conserved charges which arise as [15] 

(2n)H 



dx h n 



dx h r , 



2n)\\ 2n-i 



(2n-l)!! 

(2n)!! 
(2n- 1)!! 



In— 1 

Tro-iL^ - 



(2n-l)!! 
(2n)U 



(2n-l)!! 



2n-l 

dx tr Res L 2 



2n-l 

dx tr Res a\L 2 



dxtiU n , 

dxti o x XJ n 



(13) 



where "tr" denotes the trace over matrix indices, "Res" (Residue) stands for the coefficient of the 

p^ 1 term in the expansion and n = 1, 2, 3, (Alternatively, we can think of the two sets of charges 

as the traces of powers of L\ and L2 respectively.) Explicitly, the two sets of conserved densities 
have the closed forms 



2 M 

fc=0 x 

m 

k=l 



n 
2k 



u n-2k v 2k 



n 

2k - 1 



n n-2fc+l t ,2fc-l_ 



(14) 



The conserved quantities are, of course, defined up to overall normalization constants. Rescaling 

1 



{H n , H n 



(H n , H n ), 



n > 1, 



n(n — 1) 

we see that the equations ( |Tl|) and (|l2| ) can be written in the Hamiltonian form respectively as 



(15) 



U t 



Ut 











5U 
5U 



(16) 



The conserved densities of the system are easily seen to satisfy the usual relations for the two 
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component hydro dynamic systems |jl4| , namely, 



d 2 h n _ d 2 h n d 2 h n _ d 2 h n 

du du dv dv ' du du dv dv 

This system of equations is integrable. However, we will not get into the details of this. Our interest 
has been only to discuss it briefly so that we can bring out the similarities of this system with the 
non-linear sigma model. 



3 Lax equation for the non-linear sigma model: 

The non-linear sigma model is not exactly a dispersionless system. Consequently, we do not expect 
a Lax function in the phase space to describe such a system. We find that, like the Toda lattice, the 
non-linear sigma model does not have a scalar Lax description, but can be given a Lax description 



much like the Toda lattice [10, 13 1 in the following manner. 



Let us consider the symmetric 2x2 (multiplicative) Lax operator of the form 

/ T a lb _i_ T a jb T a lb _j_ Ta jb \ 
^ab _ / ^O^l + ^l^O \ 



S" 



where 



\ T a Tb _i_ T a lb Ta jb i ja jb j 
\ J Q J 1^~ J 1 J Q J J ~>~ J l J l / 

(Jo Jo + Ji J\) I + ( Jo a Jl + Ji Jo) vi 

L a L b , (18) 



L^Jgl+Jfa^ " " . (19) 




Let us also introduce an anti-symmetric 2x2 matrix operator 

. , . . N / 5 ab 8 + f abc J{ \ 

Tjab = f 5 abQ + f abcjc\ I ^ 

v 7 V 5 ab a + / abc Jf / 

Then, it is easy to check that the Lax equation of the form of the Toda lattice, namely, 

Qgab 



[S,U] ab = -(S ac U cb -U ac S cb ) , (21 



dt 



leads to the equations for the non- linear sigma model (|l]). (U ab , here, is an operator and should not 
be confused with the matrix U in (|j) of the last section, which describes the dynamical variables of 
the system much like S ab .) Parenthetically, we would like to remark here that the correspondence 
with the Toda lattice, however, is not quite complete as the operator U ab cannot be obtained as 
the Frechet derivative of the equations for the sigma model [[jl], |T^]. Just like the two component 
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Riemann equation of the last section, it is clear that 



[a 1 ,L a ] = 0, a u S' 



-tab 







<T1,U 



ab 







(22) 



which defines an analogous symmetry in the non-linear sigma model. 

Equation (pl| ) defines the Lax equation for the non-linear sigma model. As in the Toda lattice, 
it follows that the conserved densities of the system can be obtained as traces of powers of the Lax 



operator S |T0|, [13]. In fact, because of the symmetry (|22| ) in the system, it is trivial to check 
that there are two infinite sets of local conserved densities in the system, namely, 



hr. 



1 

in 



TiS n 



K 



l 

4n 



Trai5 r ' 



(23) 



Here, "Tr" denotes trace over the matrix indices as well as over the internal indices. Using the 
definition (|l^) as well as the cyclicity property of the trace, it is easy to see that we can write the 
conserved densities in the simpler form 



h n = ^tr(L 2 ) n 
4n 



(24) 



where 



L 2 = L a L a = (J 2 + J 2 )I + 2(J Ji)^ , 



(25) 



and "tr" denotes, as before, trace over matrix indices. In (25) the internal indices of the terms 
inside the parenthesis are summed. The similarity of these with the conserved densities of the two 
component Riemann equation in ( |i"3| ) is remarkable (we have chosen a particular normalization 
for simplicity). The first few conserved densities for the non-linear sigma model have the explicit 
forms, 



hx = \ (J 2 + Jl) , 



h2 = \(4 + Jlf + {JoJi) 2 , 

^3 = I (Jo + Jff + 2 (J 2 + Jl) (JoJif , 



hi = (JqJi) , 

~h 2 = (J 2 + J 2 ) (J Jl) 



h 3 = (J 2 + Jl) 2 Wi) + \Wif 



(26) 



We see that the local conserved densities in ( |24| ) do not involve spatial derivatives, much like the 
conserved densities of dispersionless systems. It is also clear now that the two sets of conserved 
charges can be mapped into those of the two component Riemann equation fli~3| ) or (|J) (up to 
normalization constants) with the identifications 

u = 4 + 4 , v = 2J Ji . (27) 

With this mapping, even the equations for the sigma model (Q) go over to those of the two com- 
ponent Riemann equation ( |i~2| ) for n = 0. However, such a map is not one to one and, therefore, is 



6 



not very useful. Nonetheless, it brings out an interesting relation between a dispersionless system 
and the non-linear sigma model. 

As in hydrodynamic systems, in the case of the non-linear sigma model, it can be shown that 
any local conserved density, Q (without spatial derivatives), must satisfy 



d 2 Q d 2 Q 



dJgdJg dJfdJl ' 



(28) 



which is the analogue of (|T7|), The conserved densities in ( |2"4| ) can be seen to satisfy this. However, 
from the form of the conserved densities in (|24|), it can be checked that they satisfy additional 



relations which can be written collectively as 

d 2 K d 2 h n d 2 h n d 2 h n 



dJ$ 84 d Jf dJl ' d 4 d4 dJf dJ\ ' 

d 2 K d 2 K d 2 h n d 2 h n 



3434 &4 ajf ' ajgaji 34 d jf 

ja dh n , ja dh n „, T a — + T a — - In h 

8h n dh n ~ ja^}n ja ^}n _ „ , 

f abc (4 3 ^ + J*Uo, f abc (4 d ^ c + jf K ^ - 



cabc ( rb dh n . T b dh n \ _ fab c ( .b^hn jh 3h n 



Tb^'"n T b"' u n \ n raoc jo ~ -» , to 

J °dTz +Jl dT§) =0 > f [ J °dTf +Ji d4 



(dK\ 
84 

dh n 
V dJf ) 



( dh n \ 

84 

8h n 
V dJf j 



(29) 



These relations are quite useful in proving the complete integrability of the system as well as various 
other interesting properties, as we will see. 

Once the two sets of local conserved charges have been obtained, it is straightforward to show 
that the non- linear sigma model has two Hamiltonian structures. Defining 



vf 




5 ab d 
-f 



abc jc 
J 




(30) 



it is easy to see that the non- linear sigma model equations, (|l|), can be written in the Hamiltonian 
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form as 





= vf 


{ * ) 





( ™± \ 

64 
SHi 

V J4 ) 





t 


8H X 


\ 


= vf 




54 








SHi 




) 




54 


) 



(31) 



In fact, using the properties in (29), it is easy to show that all the higher order equations in the 
hierarchy can also be written as 



I J o,t 



V 4t 



V 



ab 



54 
SHn 

V 84 ) 



V 



ab 



( 8H n \ 
54 

V 64 J 



(32) 



The anti-symmetry of the Hamiltonian structures in (30) is manifest and it can be seen, using 
the method of prolongation jig], that they satisfy Jacobi identity. Therefore, these are genuine 
Hamiltonian structures of the system. Furthermore, it is also easy to check, using the method of 
prolongation, that 



(33) 



satisfies Jacobi identity, where £ is an arbitrary constant parameter. As a result, the system 
has two compatible Hamiltonian structures and is truly a bi-Hamiltonian system. However, it is 
worth noting that, unlike in conventional integrable systems, here the two Hamiltonian structures 
act on different families of charges to yield the same dynamical equations. Since the criterion 
of Magri || applies to systems where the bi-Hamiltonian structures act on the same family of 
conserved charges, it is not a "priori clear that the integr ability of the system follows from the 
existence of a bi-Hamiltonian structure. This is a question that we will take up in the next section. 
We would like to point out here that under the action of V±, the conserved charge Hi generates 
spatial translations, as does Hi under the action of V%. This is another interesting feature of this 
system, namely, the "Hamiltonian" and the "momentum" (with respect to a given Hamiltonian 
structure) seem to be distributed into distinct families of conserved quantities in the non-linear 
sigma model. 



4 Complete integrability of the system: 

Before proving complete integrability of the non-linear sigma model, let us derive the recursion 
operators that relate the conserved charges of the two infinite sets. Using the relations in (^), it 
is straightforward to show that the Lax operator, S ab , defines a recursion operator for the system. 
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It can be explicitly checked that 



/ 5H n+1 \ 



SJg 



SH; 



n+l 



sr 



ab 



V 5J? J 



( ^ElL \ 
SJ b 

SHn 

V 8J\ J 



( SH n+l \ ( SH n \ 

= s ab 



SJg 



5H, 



n+l 



V 5J? J 



SJ b 



SHn 

V 5J\ j 



(34) 



Namely, the charges within each of the two sets are related recursively by the Lax operator, S a 
itself. 

Furthermore, let us define 



S" 



ab 



cb 



5 ab ai + d 



-1 



jabcjc jabc jc 





(35) 



Using (^]), it is easy to check that the second term on the right hand side of ( pop gives zero acting 
on 



It follows then that 



/ SHn \ 
5J b 

SHn 

\ 5Jl J 



or 



/ SH n \ 

s4 

SH n 

V sj\ ) 



(36) 



/ SH n \ 



5J$ 

SH n 

SJf 



06 



/ SHn \ 

SJ b o 

SHn 

V sj\ J 



( 5 Jh \ 



(?1 



V 



SJg 

SHn 

SJf 



(37) 



This is, in fact, already noted (the last relation) in (^). It is worth noting here that the recursion 
operator, (|35|), conventionally constructed from the two Hamiltonian structures in an integrable 
system, in this case relates charges of the two distinct sets. In fact, on the space of the gradients 



of the Hamiltonian, (36), it is easy to verify that 



I . 



(38) 



Namely, although the square of this recursion operator is not identity, its effect on the space ( p6| ) 
is that of the identity operator. This is quite important in that it implies that there are only two 
infinite sets of local conserved charges (without derivatives) associated with the system. 

Thus, we see that, in the case of the non-linear sigma model, there exist two recursion operators 
- S ab relates charges within a family whereas S ab relates charges between the two families of 
conserved charges. By construction S ab is related to the two Hamiltonian structures of the system. 
One can similarly ask if there is yet another Hamiltonian structure within a given family of charges 
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constructed in the conventional manner as 



v ab = V\ c S cb 



(39) 



We have checked that this structure does not satisfy Jacobi identity. (In fact, this structure does 
not even have the required anti-symmetry property. If it is anti-symmetrized by hand, the new 
structure does not lead to the correct dynamical equations, in addition to the fact that the Jacobi 
identity does not hold.) This, therefore, remains an open question. 

To prove complete integrability of the system, we note that 



{H n , H m } 1 — J 



il.r 



dh n ( \dh 



d- 



+ 



dh n ( dh 



d- 



eabc jc^m 



rabc jc 
J J 



(40) 



From the forms of the conserved charges in (24), it is easy to see that the term on the right hand 
side involving the structure constant does not contribute. As a result, we have 



{H n , H m } 1 — J 



dx 
dx 



dh n dh m dh n ~dh m 



dJg dJt dJf dJg 

P J 0,x + U J l,x 



(41) 



where 



pa 



G a 



dh n d 2 h m dh n d 2 h 



+ 



d 4 dJg 8J\ dJ\ dJg 84 ' 



dh n d 2 h r , 



+ 



dh n d 2 h r , 



&4 dJf 8J\ d,J\ dJf d 4 



(42) 



It is straightforward now to check, using (29), that 



84 



(43) 



This implies that there exists a function X such that 



pa 



dX 

dJg 



G" 



dX 



d 2 X 



8 2 X 



dJf ' 8Jgd4 84dJf 



(44) 



Substituting this into (pl|), we obtain 



d,T 



dX 

(It 







(45) 



with the usual assumptions on the asymptotic fall off of the fields. In a similar manner, it can be 
shown that 



(46) 
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Namely, all the charges are in involution with respect to the first Hamiltonian structure. The 
involution with respect to the second Hamiltonian structure is also easy to show along similar lines, 

{H n ,H m } 2 = = ^H n ,H m \ = ^H n ,H m \ . (47) 

This proves that the non-linear sigma model is completely integrable. 



5 Discussions: 

In this letter, we have proposed a Lax equation for the non-linear sigma model which directly 
generates all the conserved local quantities (without any derivative) of the system. There are, in 
fact, two infinite sets of conserved charges and the similarities of this system with the two component 
Riemann equation are discussed. We have obtained the two Hamiltonian structures of the system 
and have shown that they are compatible thereby proving that the system is bi-Hamiltonian. 
However, unlike conventional integrable systems, here the two Hamiltonian structures act on charges 
in the different sets to give the same dynamical equations. As a result, integr ability of the system 
does not follow in spite of the fact that it is a bi-Hamiltonian system. We have constructed two 
recursion operators for the system — one that relates conserved quantities withing any given family 
and the other that relates the charges between the two families. We have explicitly shown that the 
two sets of charges are in involution with respect to either of the Hamiltonian structures, thereby 
proving the integrability of the system. We have also brought out several other interesting features 
associated with the system. 

In the conventional discussion of the non-linear sigma model ||| , there is a standard Hamil- 
tonian structure given by 

/ f abc JR 5 ab d + f abc Jf \ 
^std =\ , \ ■ (48) 

V $ ab d + f ahc J{ ) 

This is manifestly anti-symmetric and can be easily checked to satisfy Jacobi identity. On the other 
hand, this seems to be quite different from the two Hamiltonian structures that we have obtained 
in ([}(]). We note that the Hamiltonian structure in (^) acting on H\ leads to the non- linear sigma 
model equations, (||). On the other hand, we note that 

(fabc jc fabc jc \ 
r4 A ) ' 

Using (p9|), it is easy to see that this difference gives vanishing contribution on the space (|36|). As 
a result, on the space (|36|) , the two structures T>® b d and T>f b have identical effect. (It is interesting 
to note that, in spite of this, the two structures, T> ab and T>^ b d are not compatible.) This feature 
is completely new, namely, we have two Hamiltonian structures, satisfying Jacobi identity, which 
give the same equations acting on the same Hamiltonian (normally multi-Hamiltonian structures 
act on different Hamiltonians) . It is not a priori clear, therefore, how one can select between these 
two structures. The difference in the two Hamiltonian structures lies in their action on the space 
of charges that are not gauge invariant (which have a free internal index). Namely, we know that 
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the non-linear sigma model has an infinite set of conserved non-local charges, the first of which is, 
in fact, local and is of the form 



This charge is, in fact, supposed to generate the global rotations of the system through the action 
of the Hamiltonian structure of the system. It can be checked that the action of T>^ d acting on this 
leads to the correct rotation while Df b does not. This may be a reason to choose the conventional 
Hamiltonian structure over Vf b . These are, however, open questions that need further study. 
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